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ABSTRACT 

This paper reports the spectral and timing analyses of the quiescent low-mass X-ray binary 
(qLMXB) U24 observed during five archived Chandra/ ACIS exposures of the nearby globular cluster 
NGC 6397, for a total of 350 ks. We find that the X-ray flux and the parameters of the hydrogen 
atmosphere spectral model are consistent with those previously published for this source. On short 
timescales, we find no evidence of aperiodic intensity variability, with 90% confidence upper limits 
during five observations ranging between <8.6% rms and <19% rms, in the 0.0001-0.1 Hz frequency 
range (0.5-8.0 keV); and no evidence of periodic variability, with maximum observed powers in this 
frequency range having a chance probability of occurrence from a Poisson-deviated light curve in excess 
of 10%. We also report the improved neutron star (NS) physical radius measurement, with statistical 
accuracy of the order of ~10%: i?NS = 8.9^g'gkm for Myss = 1.4 M©. Alternatively, we provide the 
confidence regions in mass-radius space as well as the best- fit projected radius i?oo = 11.9^;^;° km, as 
seen by an observer at infinity. The best-fit effective temperature, kTes = 80^5 eV, is used to esti- 
mate the neutron star core temperature which falls in the range Tcoro — (3.0 — 9.8) xlO^K, depending 
on the atmosphere model considered. This makes U24 the third most precisely measured NS radius 
among qLMXBs, after those in wCen, and in M13. 
Subject headings: stars: neutron — X-rays: binaries — globular clusters: individual (NGC 6397) 



1. INTRODUCTION 

The emission from low-mass X-ray binaries in quies- 
cence (qLMXBs) is routinely studied to provide useful 
constraints on the physical models of the interior of neu- 
tron stars (NSs). The low luminosity (10^^-10^^ erg s~\ 
4-5 orders of magnitude lower than the outburst lu- 
minosities) of these objects was first observed in the 
post-outburst stages of the transient LMXBs Cen X- 
4 and Aql X-1 (van Paradijs et al. 1987), and initially 
interpreted as a thermal blackbody emission powered 
by some low-level m ass accretion onto the NS surface 
(|Verbunt et al.lll994[ ). 

In an alternate explanation for the energy source of 
qLMXBs, the luminosity is provided by the energy re- 
leased during accretion episodes by pressure-sensitive nu- 
clear reactions (electron captures, neutron emission, and 
pycnonuclear reactions) in the NS d eep crust. The the- 
ory of deep crustal heating (DCH; I Brown et al.l Il998| ) 
describes how the accreted matter piles up at the top of 
the NS surface, forcing the matter underneath to deeper 
layers of the crust, and how the energy is deposited in 
the crust. The resulting nuclear re action chain releases 
~ 1. 5 MeV per accreted nucleo n (see lGupta et al.l ()2007l ) 
and IHaensel fc Zdunifl (J2008D for details about crustal 
heating models), and gives rise to a time-average lumi- 
nosity directly proportional to the time-averaged accre- 
tion rate: 
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where Q is the average heat deposited in the NS crust 
per accreted nucleon. 

It was also suggested that the observed thermal spec- 
trum of qLMXBs is the result of the energy deposited 
in the crust, heating the NS core during outbursts and 
thermally re-radiating away from the crust through the 
NS a tmosphere on core-cooling timescales (jBrown et al.l 

Il998f l. It is thought that this atmosphere is com- 
posed of pure hydrogen since heavy elements gravita- 
tionally se ttle on short ti mescale (~ seconds, iRomanil 
1987; Bild sten et al.|[l992l ) once the accretion from the 
low-mass evolved companion star onto the NS surface 
shuts of f after an o utburst. Mod e ls of NS H atmo- 
sphere (IRaiagopaT fc Romam] |1996l: [Zavlin et all 119961 : 

IMcClintock et al.l 12004 iHeinke et al.l 120061 ) . now rou- 
tinely used for qLMXBs, imply emission-area radii con- 
sistent with the entire area of the NS, compared to 
the ^ 1 km emission-area radii in the previously im - 
posed blackbody approximation (jRutledge et al.l [19991) . 
Spectral fitting of qLMXBs with such models also leads 
to the determination of the projected radius (as ob- 
served from infinity) defined by Roo = Rns97^ = 

Rns (1 — 2G'Mns/^nsc^) , where i?NS is the phys- 
ical radius of the NS. 

In the spectra of some qL MXBs in the field of the 
Galaxy (f or example Cen X-4.lRutledge et al.ll2001bL and 
Aql X-1, iRutledge et al.l r2001al ). an additional power- 
law component, dominating the spectrum above 2keV 
and unrelated to the H-atmosphere thermal emission, 
is observed. Proposed interpretations of this power- 
law tail include residual accretion onto the N S mag - 

netosphere (jGrindlav et al.ll2001at iCackett et al.ll2005f l. 

ill.cal 



shock emission via the e mergence of a magnetic field 
(jCampana fc Stellal 120000 . or an intrabinary shock be- 
tween the winds from the NS and its companion star 
(jCampana et al.ll2Q04t ). However, recent analyses of the 
quiescent emission of LMXBs have shown that variations 
in the non-thermal component are correlated to the vari- 
ations in the thermal component. This suggests the pres- 
ence of a variable low-le vel accretion on the NS (for the 
LMXB XTE J1701- 462.lFridriksson eTaI| [20Tol and for 
the LMXB Cen X-4, iCackett et al.ll2010f ). 

Another characteristic of qLMXBs is the expected 
lack of strong variability on long or short timescales 
since their emission is dominated by the thermal ra- 
diation from the interior of the NS. While the ther- 
mal component is not expected to display intensity 
variability, other emission mechanisms (like those men- 
tioned in the previous paragraph) may be responsi- 
ble f or intensity and spectral variations (jBrown et al.l 
|1998[ ). However, recent outburst episodes can gener- 
ate variations in the quiescent thermal luminosity on 
days to years timescales, as obser ved for the LMXB 
KS 1731-260 (iRutledEe et al."2002", following the mod- 
els described in Ushomir skv & Rutlcdgc 2001, see also 
iBrown fc Gumming! 120090 . " 

The DCH/H-atmosphere interpretation has been 
applied to a large number of historical transient LMXBs 
and provided Roo measurements from high signal-to- 
noise (S/N) spectra. The following list i s, to the best of 
our k nowledge, exhaustive: 4U 1608-522 (Rutlcdge et al.l 
[l999l), 4 U 2129 -^47 (Rutlcdgc et al. 2000) Cen X-4 
ICampan a et al.l 12000: Rutlcdge et al.l l2001bO. Agl X-1 
(JRutledge et all 12001a,) . KS 1731-260 (iRutledge et al.' 
I2002D . XTE J2123-058 ( Tomsick et al. 2004), 
EXO 1747-21 4 (iTomsick et al.l l2005l). M XB 1659-29 
(ICackett et all [2006|) , IM 1716-315 (| Jonker et all 
I2007a0. IH 1 905-^ 000 ITonker et al.ll2007bl). 2S 1803-245 
(iCornelisse eFall 120071). 4U 1730-22 (ITomsick et all 



TABLE 1 
Chandra Observations of NGC 6397 



2009D, and 



120071). EXO 07 48-676 (IDegenaar et all 
XTE J1701-462 (jFridrikss on et al.ll2010[ ). 

However, the 10%-50% systematic uncertainty in the 
distances to qLMXBs in the field directly affects the i?oo 
measurement uncertainty. Obtaining precise constraints 
on the dense matter equation of state (EoS) is the ob- 
servational motivation for measuring the radii of NSs, 
requiring ~ 5% accuracy on R^c, t o be useful for thi s 
purpose (jLattimer fc Prakashll200l iSteiner et al.ll2010l ). 
The known distances to globular clusters (GCs ) and their 
expec ted overabundances of X-ray binaries (jHut et aD 
Il992f ) has motivated the search for qLMXBs in the core 
of GCs. 

T here are 26 G C qLMXBs spectrally identified so far 
(see iHeinke et all [2003 and Guillot et al. 2009 for two 
complementary lists). However, some of them have 
poorly constrained radius and temperature measure- 
ments. This can be due to low count statistics and/or a 
large amount of galactic absorption (for example A^h = 
1.2 X 10^^ atoms cm~^ for the GC Terzan 5) in their direc- 
tion which alters the low-energy end (0.1-1 keV) of the 
spectra where the H-atmosphere spectrum of qLMXBs 
peaks. Therefore, their identification is regarded as less 
secure. Most qLMXBs in GCs have low-S/N spectra, and 
therefore, have rather large uncertainties on their Roo 
measurements (~ 15% or more), precluding their use for 



Obs. 


Starting 


Exposure 


Detector 


Mode 


ID 


Time (TT) 


Time (ks) 






79 


2000 Jul 31 15:31:33 


48.34 


ACIS-I3 (EI) 


E 


2668 


2002 May 13 19:17:40 


28.10 


ACIS-S3 (BI) 


E 


2669 


2002 May 15 18:53:27 


26.66 


ACIS-S3 (BI) 


E 


7460 


2007 Jul 16 06:21:36 


149.61 


ACIS-S3 (BI) 


VF 


7461 


2007 Jun 22 21:44:15 


87.87 


ACIS-S3 (BI) 


VF 



Note. — All observations were performed with a 3.24104 s frame rate. 
TT refers to terrestrial time. FI and BI refer to front-illuminated and 
back-illuminated detectors, respectively. 

EoS constr aints. So far, only a few qLMXBs (in a; Cen, 
IGendre et a l. 2003 al in M13. IGendre et al.l l2003bl and 
X7 in 47 Tuc, Hcin ke et al.l l200g) have spectra with S/N 
high enough to provide useful constraints on the den se 
matter EoS (|Webb fc Barretil2007t ISteiner et al.ll20Tot) . 

The close proxim ity of t h e glob u lar cluster NGC 6397 
(d » 2.5 kpc. I jHarrisI [1991 IHansen et all [20071 : 
iStrickler et al.l 120091 ) , the moderately low galactic ab- 
sorption-^ (A^H = 0.14x10^^ atoms cm^^, noted A^h,22 — 
0.14 hereafter) in its direction makes it a useful target 
for the spectral identification of qLMXBs. The discov- 
ery observation of the qLMXB U24 ^ in NGC 6397 pro- 
vided modest constraints on the NS projected radius: 
R^ = 4.9^1"* km (Grindlav et al.ll200lU GrOl hereafter). 
The proximity of U24 to the GC core requires the use 
of the Chandra X-ray Observatory's angular resolution 
to positionally and spectrally separate the qLMXB from 
other sources in the crowded GC core. U24 lies at 
dc — 6.8 Tc w 20" from the GC center (core radius 
re = 0.05' and half-mass radius mivi = 2.33', NGC 6397 
is a core-collapse cluster). The reported effective tem- 
perature of U24 was fcToff' = 57-92 eV (90% confidence 
interval, GrOl). 

This paper presents the spectral and timing analyses 
of five archived deep Chandra- AGIS observations of U24, 
located in GC NGC 6397. These long exposures provide 
the high-S/N data necessary to confirm the qLMXB clas- 
sification of U24 by obtaining precise i?NS estimation. 
The lack of variability also supports this classification. 
The organization of this paper is as follows. Section [5] 
describes the data reduction and the analyses. Section [3] 
presents the results of the analyses. Section |4| provides a 
discussion of the results and Section [5] is a short conclu- 



2. DATA REDUCTION AND ANALYSIS 

2.1. Observations, Source Detection, and Count 
Extraction 

We analyze one archived Chandra / ACIS-I and four 
archived Chandra / ACIS-S observations of NGC 6397 
(Table [Ij. The source detection and the data anal- 
ysis are performed using the CIAO v4.1.1 package 
(jFruscione et al.l[2006[ ). The event files (level=l) are re- 
proce ssed with the latest calibration files from CALDB 
v4.1 (jGraessle et al.l [20071 with the latest effective area 



^ From http://cxc. harvard.edu/toolkit/colden.jsp using the 
NRAO data (iDickev &i Lockman 199lJ}. 

■^ The source name U24 was used in the discovery paper (GrOl) 
and will be used throughout this work. 
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TABLE 2 
Position of U24 in Chandra Observations of NGC 6397 



Obs. ID 



R.A. 
(J2000) 



{") 



Decl. 
(J2000) 



(") 



Detection 
Significance 



Reference 



79 


17h4om41.421= 


0.02 


-53°40'04.73 


0.02 


- 


GrOl 


79 


171140^41.4595 


0.6 


-53°40'04.47 


0.6 


128(T 


This work 


2668 


17^40^41.489= 


0.6 


-53°40' 04.38 


0.6 


152o- 


This work 


2669 


17^40^41.485= 


0.6 


-53°40'04.53 


0.6 


149(T 


This work 


7460 


17h40™41.486= 


0.6 


-53°40'04.60 


0.6 


302(T 


This work 


7461 


17h40™41.488= 


0.6 


-53°40'04.54 


0.6 


228(T 


This work 



Note. — Ar^ a. and AdccI. arc the uncertainties on the position, dominated by the Chan- 
dra sy stematic uncertainty in this work. The positions reported previously (Grindlay ct al.l 
I2001bl , noted GrOl) have been corrected for the Chandra systematic uncertainty (sec details 
in Section l3.H . 



maps, quantum efficiency maps, and gain maps), as rec- 
ommended in the CIAO Analysis Thread ^^Reprocessing 
Data to Create a New Level=2 Event File" to include the 
up-to-date CTI corrections (charge-transfer inefficiency). 

The re-processed event files are analyzed including 
counts in the 0.5-8.0 keV energy range. The full-chip 
light curves do not show evidence of background flares in 
any of the five observations, allowing use of the full expo- 
sure time of each observation. The low flux of the source 
of interest, U24, allows neglecting pile-up; the count rate 
of ~ 0.06 photons per frame corresponds to a pile-up 
fraction of less than 2% •^. 

For each observation (ObsID), the source detection is 
performed with the CIAO wavdetect algorithm. An ex- 
posure map is created using the task mkexpmap prior to 
the source detection. The wavdetect exposure threshold 
expthresh is set to 0.1 and the wavelet detection scales 
are set to scales = "1.0 2.0 4-0 8.0". Thirty-five sources 
are detected (cr > 3) on the ObsID 2668, 37 sources on 
ObsID 2669, 66 sources on ObsID 7460, 48 on ObsID 
7461, and 37 on the ACIS-I observation ObsID 79. 

This paper is solely focused on radius measurement 
and timing analysis on the qLMXB U24. While the 
source detection is performed over the whole ACIS chip, 
the following analysis pertains only to the source U24. 

Counts are extracted with the CIAO script psextract 
around the source position in a circular region of ra- 
dius 3", which ensures that 98% of the enclosed en- 
ergy fraction at 1 keV is included ^. The closest de- 
tected source, located at 10.6" distance from U24, has 
21.6 counts (background subtracted) within 1.5". It con- 
tributes to <^ 0.04 contamination counts within the ex- 
traction radius of U24 (on the longest observation). The 
background is extracted from an annulus centered at the 
qLMXB position with inner radius 5", and outer radius 
30". Other detected sources within the background an- 
nulus are excluded with a 5" radius region, which elimi- 
nate 99.8% of source counts in the background region. 
For the deepest observation (ObsID 7460), 15 counts 
from other sources are within the extracted background 
(which contains 6187 counts). In other words, these con- 
straints ensure that ~ 0.25% of the background counts 
are due to other sources. Finally, the extraction radius 

^ Chandra Observatory Proposer Guide vl2.0, Figure 6.18, De- 
cember 2009 

** Chandra Observatory Proposer Guide vl2.0, Figure 6.7, De- 
cember 2009 



(containing 98% of the ECF) does not require to apply 
a correction to the fiux. 

Following the CIAO Science Thread "Creating ACIS 
RMFs with mkacisrmf", the response matrices files 
(RMFs) are recalculated prior to the spectral analysis 
since the RMFs obtained from psextract are not suited 
for ACIS observations with focal plane temperature of - 
120°C (the usual mkrmf command does not use the latest 
calibration available in the case of -120°C ACIS imaging 
data). In addition, the ancillary response files (ARFs) 
are also recalculated using the energy grid of the newly 
obtained RMFs. Overall, the extracted spectra, together 
with the RMFs and ARFs, are used for the spectral anal- 
ysis. In those spectra, the effect of background counts 
can be ignored. Indeed, in the worst case (for ObsID 
7460), the number of expected background events ac- 
counts for 2.4% of the total number of counts in the ex- 
tracted region (78.0 background counts out of a total of 
3188 counts), so that the background is neglected for the 
spectral analysis. 

2.2. Spectral Analysis 

For each of the five observations, two spectral files are 
created, one with unbinned events (for fitting with the 
Cash-statistics, Cash ,1971) and one with binning (for 
fitting with the x^-statistics). For the latter, the bin 
width in the 0.5-1.5 keV energy range matches the energy 
resolution of the ACIS-S3 chip, i.e., ~0.15keV. Above 
1.5 keV, four wider bins (0.3 keV, 0.6 keV, and two 3keV 
wide spectral bins) are created. In some cases of low 
count statistics, the last 2 or 3 bins are grouped together 
to maintain a minimum of 20 counts per bin. The main 
criterion for the creation of the spectral bins is the energy 
resolution of the detector, but the 20 counts minimum 
is imposed to ensure approximate Gaussian uncertainty 
in each bin. Such a binning avoids an artificially small 
reduced-x^, and conserves the validity of x^-statistics. 

Spect ral fitting is pe rformed with the software XSPEC 
vl2.5.1 ()Arnaudlll996[ ) using the publicly available model 
of NS H-atniosphe re nsatmos ()McClintock et al.l [200l 
'Heink e et al.l 120061 ) . The model assumes non- magnetic 
NSs and has been computed for a range of surface grav- 
ity g — (0.1 — 10)xlO^''cms~^. For the normalization 
parameter, nsatmos uses the emitting fraction of the NS 
surface. It is kept fixed to unity in this work; in other 
words, the whole NS surface emits. The distance pa- 
rameter is held fixed as well at the value of NGC 6397, 



d = 2.5 kpc ([Hansen et al.ll2007t iStrickler et al.ll2009[ ). 
The NS mass is assumed to be 1.4 Mq. Finally, the galac- 
tic absorption is taken into account using the phabs mul- 
tiplicative model, with iVfj, the hydrogen column density 
parameter, set to -/Vh.22 = 0.14. The errors on the best- 
fit parameters (-RnSj fcT'off) are calculated using the com- 
mand error in XSPEC with 90% confidence or using the 
command steppar. Confidence contours in mass-radius 
space are obtained with the steppar command with both 
the mass and the radius as free parameters. The results 
of the spectral analysis are presented in Section 13.21 

2.3. Variability Analysis 

For each of the five observations, we perform two anal- 
yses to search for source variability on timescales shorter 
than the duration of the observations, and one analysis 
for timescales spanning the time between the first and 
last observations. 

Power Density Spectrum (PDS). The data from 
each observation were binned into a discrete light curve, 
with time bin size equal to the time resolution used in 
the observation (AT = 3.24104 s). A standard fast- 
Fourier transform (FFT) of t he discrete data into fre- 
quency space was produced fPress et al. 1995), using 
an open-source FFT algorithm ( Frigo fc Johnson 1998 ). 
This produced a Fourier transform of the data, covering 
the frequency range 1/T — 0. 15427Hz, where T is the du- 
ration of the observation, and 0.15427 Hz is the Nyquist 
frequency, with discrete frequency resolution 1/T, and a 
total of N/2 frequency bins, where N = T/AT is the 
number of time bins in the light curve. The resulting 
Fourier transform was then converted into a PDS, where 
the power Pj in each frequency bin j is Pj = a^ + b'j, 
where Uj and bj are the real and imaginary parts of the 
Fourier component associated with a frequency fj = j /T, 
producing a PDS. The data were then norma hzed accord- 
ing t o a standard prescription for analyses (jLeahv et al.l 
fl98l . 

Short-Timescale (<lday) Variability. Short- 
term variability in each ObsID is assessed by visual in- 
spection of the five sources light curves . In add ition, 
a Kolmogorov-Smirnov (K-S) test (Press et al.lll995i) is 
performed to quantify the consistency of the temporal 
distribution of counts in each observation with a constant 
count rate. More specifically, the integrated-count light 
curve is compared to a linear distribution using a K-S 
test. A low K-S probability would indicate that the inte- 
grated light curve is significantly different from a linear 
distribution, demonstrating the presence of variability on 
the timescale of the observation. 

Long-Timescale (^^months-years) flux variabil- 
ity. Long-term variability is investigated to determine if 
the flux remained constant over the course of the five ob- 
servations, between 2000 and 2007. This is performed by 
adding a multiplicative constant parameter to the spec- 
tral model, and by fitting the five spectra simultaneously. 
The constant is kept fixed at a value c = 1 for one spec- 
trum (ObsID 79) and as a free and untied parameter 
for the remaining four spectra while the remaining spec- 
tral parameters (iVn, -RnSj kT^s) are assumed to be the 
same across all observations. Best-fit c values statisti- 
cally consistent with unity would demonstrate that the 
source flux remained constant on the timescale of the flve 
observations. 
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Fig. 1. — Folded spectra of all five Chandra- ACIS observations 
of U24 in NGC 6397, in the 0.5-8.0 keV energy range. Each color 
corresponds to the spectra of each observations: ObsID 79 in black, 
ObsID 2668 in red, ObsID 2669 in green, ObsID 7460 in blue, Ob- 
sID 7461 in cyan blue [See the electronic edition of the paper for 
a color version of this figure]. The solid lines are the best-fitting 
model of an NS H-atmosphere model: NSATMQS, with A'h,22 = 0.14 
and AfjMg = 1.4 Mq kept frozen. The spectral binning was per- 
formed in accordance with the Chandra energy resolution (see Sec- 
tion [2T2ll. The fit is statistically acceptable with x?/dof (prob.) = 
1.49/45 (0.02). The residuals, shown in the lower part, indicate 
evidence of a count defic iency between 0.8 keV and 1.2 keV (see 
Figure [2] and Section 13.211 , as well as an apparent excess of counts 
at large energy, above 2.5 keV, which is fit with a power-law com- 
ponent with a = 1 to estimate the contribution to the total flux 
(see Section l3.2l l. 

3. RESULTS 
3.1. Positional Analysis 

The position of the source reported in the discov- 
ery observation (ObsID 79) is R.A.=17''40™41.421" and 
deck = -53°40'04.73" (GrOl). The authors corrected 
for the systematic Chandra uncertainty in the point- 
ing (0.6") by cross-identifying cataclysmic variables on 
Chandra and Hubble Space Telescope (HST) observa- 
tions. The positional uncertainty on the source is sub- 
arcseconds and consists in the residual error from the 
correction and the statistical uncertainty from the detec- 
tion algorithm. 

For comparison, the source detection in all five observa- 
tions presented here led to the positions listed in Table [51 
It is concluded that the position of U24 is consistent with 
that reported previously, within the 0.6" Chandra sys- 
tematic uncertainty. No correction for the systematic 
uncertainty is performed here since the identification of 
U24 is free of source confusion and since the purpose of 
this paper focuses on the spectral analysis. 

3.2. Spectral Analysis 

The spectral analysis is performed using the Cash- 
statistics on the unbinned data (neglecting the back- 
ground) and using x^-statistics on the grouped data. Re- 
sults are presented in Table [31 listing the best-fit parame- 
ters for Cash and x^-statistics, freezing and thawing A^h- 
All reported errors are 90% confidence. The x^-statistics 
null hypothesis probability confirms the viability of the 
fitted model while the Cash-statistics assumes that the 
model describes the data. Nonetheless, the use of Cash- 
statistics provides smaller uncertainties on the best-fit 
parameters, specifically for the NS radius measurements. 

A simultaneous fit is also performed using all five spec- 
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TABLE 3 

Spectral Parameters of U24 in Chandra Observations of NGC 6397 



Obs. ID 


A^H,22 


Ans (km) 


kna (eV) 


Fx 


Statistics 


PL Contrib. 


79 


(0.14) 


lO.Tiy 


731115 


1 t;n+0.14 
l-OU_Q22 


Cash; 256.7 (100%) 


- 




0.09lO;0| 


6.8l?-i 


9711 


1.27 


Cash: 254.2 (0%) 


- 




(0.14) 


lo.ill.t 


75l« 


1 49+0.11 
l-^y-0.22 


xS/dof (prob.) = 0.76/7 (0.62) 


< 8.8% 




r. ,f. + 0.05 

'J-^U~0.03 


7 9+3.2 


93lg 


1.30 


xg/dof (prob.) = 0.60/6 (0.73) 


- 


2668 


(0.14) 


7 9+^.4 
'•^-1.8 


95lfg 


1 47+^0 
^•4'_1.34 


Cash: 213.3 (100%) 


- 




13+0-03 


7 1+2.8 
'■1-1.6 


96+26 


1.44 


Cash: 213.1 (0%) 


- 




(0.14) 


7 S + S-"! 
'•^-2.8 


88l« 


-I /I4+1.20 


xS/dof (prob.) = 1.81/9 (0.61) 


< 3.9% 




13+004 


7 2+4.6 
'■^-2.2 


94l« 


1.38 


xl/Aoi (prob.) = 1.97/8 (0.46) 


- 


2669 


(0.14) 


10 0+^" 


76+^' 


-I c9+0.0b 
1-^^-0.26 


Cash: 159.0 (98%) 


- 




n+0-05 

'J-^^-0.02 


7p,+5.4 


96lS 


1.36 


Cash: 157.7 (0%) 


- 




(0.14) 


a Q+3.0 
y-y-4.9p 


77l?J 


-1 c-i+0.05 
l-''l-0.40 


xS/dof (prob.) = 0.63/8 (0.74) 


< 4.2% 




12+0-07 


7 9+8.8 
'■^-2.2p 


93l^t 


1.37 


X^/dof (prob.) = 0.66/7 (0.70) 


- 


7460 


(0.14) 


9.9±i:f 


75l^ 


1 qc + U.04 

1.00_g Qg 


Cash: 367.0 (100%) 


- 




11+0-01 


6.9lS;^ 


94l« 


1.18 


Cash: 355.1 (0%) 


- 




(0.14) 


lO.ltli 


741^ 


1-351^:^7 


xS/dof (prob.) = 2.32/8 (0.02) 


< 4.7% 




11+0-02 


7.011:^ 


931^? 


1.19 


X?/dof (prob.) = 2.14/7 (0.04) 


- 


7461 


(0.14) 


6.7^,-0 


ioo+i« 


l-'^'-1.19 


Cash: 291.0 (0%) 


- 




14+0-02 


6.6+11 


1011^4 


1.38 


Cash: 291.0 (0%) 


- 




(0.14) 


e.ilg:^ 


1041- 


1 07+1.33 

J-3'-1.07 


X^/dof (prob.) = 0.97/7 (0.45) 


< 3.3% 




n 1 7+0.06 

'J-l'-0.04 


10-7tr,, 


741?^ 


1.57 


X^/dof (prob.) = 1.03/6 (0.41) 


- 


All 5 


(0.14) 


8.9±r6 


80+^ 


1 -30+0.02 
^••^^-0.06 


Cash: 1289.9 (100%) 


- 




12+0-01 


6.91J;? 


95l?o 


1.28 


Cash: 1276.1 (69.2%) 


- 




(0.14) 


9.311:0 


78l^ 


1 oQ+0.04 
l-''^-0.06 


X^/dof (prob.) = 1.49/45 (0.02) 


< 3.7% 




12+002 


7.2l- 


921- 


1.29 


X?/dof (prob.) = 1.16/44 (0.03) 


— 


Note. — 


The model ut 


ed for all spectre 


1 fits presented 


here is the NS 


H-atmosphere model WSATHOS, with the dista 


nee to the source 


fixed at the 1 


/aluc d = 2.5 kpc and the ma 


,B of the NS fixed at Mj^g = 


1.4 M"q. Wjj 22 '^ *^® galactic absorptio 


, Nh in units of 


10^^ atoms cir 


~2. Values 1 


n parentheses ar 


e held fixed for 


the spectral fits and quoted uncertainty are 90% confid 


nee . Fj^ is the 


unabsorbcd fl 


ux in units of 


10-13 erg cm-2 


s~^ in the 0.5 


— 10 kcV range. 


The uncertainties on the fluxes were only 


calculated for the 


fixcd-JVj^ fits 


For ObalD 7460, the x^ value 


greater than 2 


prevents XSPEC to directly estimate the error region. It w 


as obtained using 


the command 


steppar. "PL c 


ontrib." refers t 


3 the upper lim 


t of a powcr-la 


w component contribution to the total flux 


(see Section|3.2|. 


The percent 


n parenthesis 


following the Ca 


sh~statistic vali 


le indicates the 


goodncss-of-fit. The best fit parameters 


and the statistics 


information p 


rovided do not 


include the pow 


er-law compone 


rt. 







tra, therefore increasing the count statistics and provid- 
ing better constraints on the best-fit parameters (Fig- 
ure [1]). Prior to that, the source long-term variabihty is 
inspected, as described in Section [^3 The best-fit val- 
ues (with Cash-statistic) for the multipHcative constants 
described in Section [Z3l are c = 1.02 



-,+0.09 



4+0.09 



0.921°;°^, and c 



n qc;+0.07 



r — 1 V\V 
0.10' ^ ~ -"^-"^-O.IO' 



c = U.y^_!_5'57, and. c = U.yt)_!_5'5g, for the observations 
ObslD 2668, ObslD 2669, ObslD 7460, and ObslD 7461, 
respectively, and c = 1 (fixed) for ObslD 79. All best- 
fit values are statistically consistent with unity (within 
1.5(t), indicating that the source flux did not vary on 
long-term timescale. This allows for simultaneous spec- 
tral fitting with the constant multiplier fixed at the value 
c = 1 for all five spectra. 

The upper limit of a power-law contribution to the 
total flux is also estimated. To do so, a power-law com- 
ponent with fixed photon index a = 1 is added to the NS 
atmosphere model and the flux of this component using 
the upper limit of the power-law normalization parame- 
ter is measured. It is found that the power-law contribu- 
tion accounts for < 3.7% of the total flux (90% confidence 
upper limit), when estimated from the simultaneous fits. 
Upper limits on the contribution of a power-law com- 
ponent for the individual spectra are also indicated in 
Table El 

There is some evidence for a count deficiency between 



0.8 and 1.2 kcV on ObslD 7460 given the assumed spec- 
trum and instrument calibration (see Figure [2]). We 
parameterize this apparent dip in the spectrum with 
a notch component. The best-fit notch central en- 
ergy is E ^ 0.96lJl:l!^ keV, and the best-fit width is 
W — 42 zL 17 eV. The added component improves the 
statistics, x^/dof (prob.) — 1.35/6 (0.23), without alter- 
ing the best fit NS temperature and radius. 

We also investigate the statistical significance of this 
deficie ncy. Using the metho d described in a previous 
work (jRutledge fc Sakdl2003f) . we estimate the proba- 
bility of observing such a deviation from a continuum 
model. First, the observed spectrum is convolved with 
the energy redistribution of the ACIS-S detector. Then, 
using a Monte Carlo approach, it is shown that the con- 
volved spectrum does not exceed the 99% confidence lim- 
its envelope obtained from 10,000 simulated spectra (Fig- 
ure [3]). They were created using the best-fit nsatmos 
model (i?NS = 8.9 km, kT^s = 80 eV, Mns = 1-4, 
d = 2.5 kpc, and Nu,22 = 0.14, see Table [S]). The maxi- 
mum deviation corresponds to a 98.2% confidence, which 
does not constitute sufficient evidence to claim the de- 
tection of an absorption line, but which should be inves- 
tigated in more details with higher S/N observations. 

For completeness, the results of spectral fits with other 
models are provided. Using an absorbed simple black- 
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body model, the fit to the spectra (aU five binned spec- 
tra) is not statistically acceptable: xj/dof (prob.) — 
3.51/45 (2x10^1'') with A^h,22=0.14 fixed, and x^/dof 
(prob.) = 1.631/44 (5x10^^) with A^h,22 allowed to vary 
A thermal bremsstrahlung model is also fit to the spectra, 
leading to a statistically acceptable fit (x^/dof (prob.) = 
1.45/45 (0.03)) with best-fit parameter kT = 391±10eV. 



Observed Convolved Speclrum - 
99% upper limit 
99% lower limit - 



1 : 

- +± \ - 



Fig. 2. — Folded spectrum of the ACIS observation 7460. As for 
Figure [l] the solid line shows the best-fitting NSATMOS model. The 
count deficiency observed between 0.8 keV and 1.2 keV does not 
exceed the 99% limit from the MC simulations of 10,000 spectra. 
In other words, the detection is not signific ant e nough to claim the 
presence of an absorption line (see Section l3.2l l. 



3.3. Variability Analyses 

We find no evidence of broadband variability as a 
function of frequency in the PDS. For broadband vari- 
ability, uncertainties were assigned to each frequency 
bin equal to the square root of the theoretical vari- 
ance in the power, approp riate to the assumed normal- 
ization (|Leahv et al.lll983l) . and derived 90% confidence 
(2cr) upper-limits on the root-mean-square (rms) vari- 
ability. The data were then rebinned logarithmically 
and the resulting PDS were fit with a model of a con- 
stant power ^ , plus a power-law component in which the 
power scales P., oc /~", and the power-law slope was 
held fixed at a = 1. This model was fit to the data us- 
ing a Leve nberg Marquardt x^ minimization technique 
([Press et al . 1995), to find the best- fit model param- 
eters. The resulting 90% confidence upper limits on 
the broadband variability, in a frequency range 0.0001- 
0.10 Hz (used for each observation, to ease comparison 
of limits, although the longer observations are sensitive 
to variability at frequencies below this range) and across 
the full Chandra/ ACIS energy range (0.5-8.0 keV) were 
: < 12% (ObsID 79), xj=0.63 (26 dof); < 19% (ObsID 
2668), x?=l-5 (27 dof);< 11% (ObsID 2669), x^=l-32 
(26 dof); <6.4% (ObsID 7460), x^=0.98 (15 dof); <8.6% 



(ObsID 7461), xS=0.98 (9 dof). 



^ The constant accounts for the expected Poisson noise power 
in the PDS, which is approximately 2 in the PDS normalization 
used; however, we observed the Poisson level for each of the PDS 
to be suppressed, to a value of ~ 1.9, due to instrumental dead- 
time effects. To correct for this, on average, in the rms variability 
upper limits derived here, the measured rms limits were increased 
by a factor (j4)/2.0, where (A) is the best-fit Poisson power level. 




Fig. 3. — Results of the MC simulations analyzing the count 
deficiency observed in Figure [2] in the energy range 0.8— 1.2keV. 
The observed spectrum (solid line) convolved with the ACIS energy 
redistribution, C{E), is shown within the 99% upper- and lower- 
limit envelopes (dashed and dotted lines, resp ectively) obtained 
from Monte Carlo simulations (see Section l3.2| l. 

We find no evidence of periodic variability. For Ob- 
sIDs 79, 2668, 2669, 7460, and 7461 respectively, we find 
maximum (normalized) powers of P,nax =21.8, 17.0, 15.3, 
17.7, and 23.7 which, with a number of frequency bins 
of 7526, 4383, 4154, 23072, and 13812, correspond to re- 
spective probabilities of chance occurrence in all cases of 
>10%. 

Visual inspection of the light curves did not reveal any 
variation in the source count rates. Short-term variabil- 
ity was also investigated in a more quantitative way by 
comparing the integrated ligthcurve with a linear dis- 
tribution. None of the five observations showed an in- 
tegrated light curve significantly different from a linear 
distribution. More specifically, the calculated K-S prob- 
abilities were 65% (ObsID 79), 62% (ObsID 2668), 95% 
(ObsID 2669), 51% (ObsID 7460) and 16% (ObsID 7461). 
Therefore, we find no evidence of intensity variability on 
the timescale of the integration times, i.e., ^ 1 day. 

4. DISCUSSION 

4.1. i?oo Calculation 

Producing realistic constraints on the dEoS requires 
obtaining values of i?oo relying on as few assumptions 
as possible. Keeping the mass fixed for the model fit- 
ting is therefore not appropriate for that purpose. We 
estimate the value of Roo by permitting A/ns to vary, 
and calculating the contours of constant model probabil- 
ity resulting from the fits in a mass-radius space (Fig- 
ure [5]). This is done using the steppar command in 
XSPEC. The 90% confidence regions of i?NS and Mns 
are obtained from the 90%-contour in the mass-radius 



space: i?NS 



9.7l|^Jkm and A/ns 



1.13 



-1-0.47 
-0.32 



Ma 



The best-fit value of the projected radius is therefore 
Roo = 11.9 km. The calculation of the uncertainties is 
complicated by the fact that the distribution of i?NS and 
A/ns is not symmetric around the best-fit values (i.e., 
not Gaussian). Moreover, i?NS and A/ns are highly cor- 
related, as shown by the crescent shape of the contour in 
M-R space (Figure |4|). Therefore, the calculation of the 
uncertainties on Roa using Gaussian normal error prop- 
agation is not valid. We describe two methods to obtain 
the uncertainty on i?oo- 
The projected radius and its uncertainties can be ob- 
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Fig. 4. — Contour plot in mass-radius space resulting from the 
simultaneous fit of the five data sets with the NSATMOS model. The 
67%, 90%, and 99% contours are shown, and the cross indicates 
the best fit (when both the mass and the radius are free to vary). 
For an NS mass of 1.4 Mq, the range of iJNg corresponds to that 
obtained and quoted in Scction l3.2l 

tained from a tabulated version of the nsa spectral 
model (jZavlin et al.lll996l ). However, this model (or nsa) 
is less adapted than the nsatmos or nsagrav models 
(|Webb fc Barreti [20071 ) because it was calculated for a 
single value of the surface gravity g = 2.43 x 10^^ cm s^^ 
while the other two models consider a range of val- 
ues. Nevertheless, the best-fit Roo value with this model 
is: Roo — 12.lj]|!)gkm (consistent with the value cal- 
culated in the previous paragraph), for a temperature 
fcToff = 76^3 eV (x^/dof (prob.) = 1.54/44 (0.02)). 

A second method to estimate the uncertainties in- 
volves geometric construction, by reading graphically the 
error region of Roo on the M-R contours (Figure S]). 
For that, we choose to use the line of constant sur- 
face gravity (i.e., constant M [R)) that goes through 
the point {R,M) — (0 km, OAf©) and the point of best 
fit in M-R space. This line intersects the 90% con- 
tour at the points iR,M) = (9.048 km, 1.045 Mq) and 
{R,M) = (10.39km, 1.245 Mg). These two points corre- 
spond to the values Roo = 11.15 km and Roo = 12.92 km 
which are, respectively, estimates of the lower and upper 
90% confidence uncertainties on Roo, assuming a con- 
stant value of the surface gravity. 

Therefore, the projected radius and its estimated 90% 
confidence uncertainties are: Roo = 11.9lQgkm. With 
the achieved uncertainty, U24 becomes the third best 
radius measurement of a NS among the population of GC 
qLMXBs, af ter the ones in ojCen (jGendre et al.|[2003af) 
and in M13 ()Gendre et al.ll2003bn . 



4.2. Error budget 

The high S/N spectra and the precise radius measure- 
ments obtained in the work presented here can be used to 
constrain the EoS of dense matter. A high precision on 
the NS radius is mandatory to exclude some of the exist- 
ing nuclear dense matter EoSs and provide the necessary 
constraints to understand the behavior of such matter. 
However, other sources of error come into play in this 
type of measurements. To quantify the total uncertainty 
on the radius measurement presented here, we estimate 
the contribution of each source of error into an error bud- 
get, including the distance to the GC NGC 6397, un- 



certainties intrinsic to the model used, systematic and 
statistical uncertainties. In those reference s where these 
uncer tainties are discussed (for example, iHeinke et al.l 
|2006[) , only two of the three uncertainties we discuss here 
(distance and detector systematics) are addressed. No 
work that we can find in the literature discusses the im- 
pact of the uncertainty in the spectral model on derived 
model parameters; therefore, we do so here. 

• The distance to the GC was recently measured us- 
ing two independent methods. The analysis of the 
CO white dwarf (WD) sequence from deep observa- 
tions in an outer field of NGC 6397 led to a distance 
of 2.54±0.07kpc (Ha nsen et al.ll2007D . More re- 
cently, using CO WDs in central regions of the clus- 
ter, the distance was calculated to be 2.34±0.13 kpc 
(jStrickler et al.|[2009l) . The weighted mean of these 
two measurements is 2.50±0.06kpc, correspond- 
ing to a 2.4% uncertainty. The unknown line-of- 
sight position of U24 within NGC 6397 accounts 
for < 0.1% of the distance uncertainty, which can 
be neglected compared to the GC distance uncer- 
tainty. 

• The calculation of spectral model NSATMOS also 
contributes to the total uncertainty on the mea- 
sured radius. However, the previous works de- 
scribing the model do not provide a discussion 
on the fractional uncertainties in intensity due to 
conver gence during the calculation of the sp ectral 
model (jMcChntock et al.ll200l IHeinke et arlr2 006') . 
Therefore, it is not possible to evaluate the errors of 
the resulting spectra. The cite d reference for sim- 
ilar models (NSA and NSAGRAV. [Zavlin et al.|[l996D 
only provides information on the temperature cal- 
culation convergence, which does not permit an es- 
timation of the uncertainty error in the modeled 
intensity as a function of energy. 

• The statistical uncertainties are those quoted in 
Table [3] (90% confidence). This includes the 3% 
systematic uncertainty of the detector calibration, 
taken into account using the "systematic 0.03" 
command in XSPEC. 

Consequently, the distance uncertainty (2.4%) is the 
only quantifiable error not taken into account in the ra- 
dius measurement obtained from spectral fitting. It is 
therefore added in quadrature to the systematic and sta- 
tistical uncertainties to obtain the total quantifiable un- 
certainty in the radius measurement. For example, the 
upper bound uncertainty limit of i?NS was 10.1% and 
is 10.4% when accounting for the distance uncertainty. 
The lower bound uncertainty limit was 6.7% and becomes 
7.1%. Consequently, the physical radius is i?NS = 8.9^0.6 
(for AfNS = 1.4 Mq) while the estimated radiation ra- 
dius is Roo = ll-9logkm, when considering the sources 
of uncertainty listed above. In conclusion, in NGC 6397, 
the distance uncertainty of NGC 6397 alone minimally 
affects the current uncertainty on the radius measure- 
ment. 

4.3. Core temperature calculation 



The best-fit temperature and physical radius of the NS 
can be used to determine the interior temperature. This 
calculation is model dependent and due to the uncer- 
tainties in the deep atmosphere composition (the depth 
of the H/He transition in particular), two different mod- 
els are considered here. The first one assumes a layer of 
helium down to a column depth y — IxlO^gcm"^ with 
a pure layer of iron underneath. The second model con- 
siders a thin layer of He down to y = 1 x 10* g cm~^ with 
a mixture dominated by rp-processes. These two alter- 
natives take into account the extremal values for the core 
temperature, for a fixed effective temperature. 

The calcul ation, described in a previous work 
(jBrown et al.l [2002) , was performed with a fixed mass 
A^NS = 1.4 M0 and a fixed radius i?NS = 8.9 km, which 
corresponds to the best-fit value. The effective temper- 
ature used was kTcs = SOtg eV. The calculation is per- 
formed down to a column depth of IxlO^^'^gcm"^, since 
the temperature is nearly isothermal in deeper layers. 

For the first model, the resulting interior tempera- 
ture (at y = 10"gcm-2) is Tcoro = (S.STtoi?) xlO^K. 
The second model leads to the value of interior tem- 
perature Tcore = (S.QStJSJg) xlO^K. Overall, if it is 
assumed that the H/He transition depth is unknown, 
the core temperature is in the range of extreme values: 
Tcore = (3.0-9.8)xlO^K. 

5. CONCLUSION 

We have performed the spectral analysis of five 
archived Chandra observations of the qLMXBs in the GC 
NGC 6397. The ^ 350 ks of integration time available 
permitted to obtain high S/N spectra and improve the 
radius measurement. More specifically, the simultaneous 
spectral fitting of all five observations, using an NS H at- 
mosphere, allowed us to provide constraints on the NS ra- 
dius with ^^10% statistical uncertainty (90% confidence). 
This confirmed the qLMXB nature of the X-ray source. 
Therefore, the measured NS properties are i?NS = 
8.9l[]|km and kT^s = 80^5 eV, for A/ns = 1.4 Mq, and 
assuming an NS with an atmosphere composed of pure 
hydrogen. The estimated interior temperature lies in 
the range Tcoic = (3.0 - 9.8) xlO^K. In the 0.5-10 keV 
range, the flux corresponding to these best-fit parameters 



is Fx = (l-39;oo6) X 10^1^ erg cm^^ g-i^ equivalent to 
a luminosity of Lx = (1.04^0.05) x 10^^ erg s"^ at a dis- 
tance of 2.5 kpc. The spectra did not show evidence for 
a power-law component as inferred by the upper limit on 
the contribution to the total flux of 3.7%. 

The results of this analysis were consistent with those 
of the discovery observation (GrOl); the reported NS 
radius and temperature were i?oo = 4.9_]^ km and 
fcTcff = 57-92 eV. No optical counterpart was detected 
on the HST observations, with a limiting magnitude of 
Mv > 11 (GrOl). 

No variability was observed on long timescales. There- 
fore, unless an outburst (for which there is no evidence) 
happened between the observations-between 2000 and 
2002, or between 2002 and 2007-we conclude that the 
source remained in its quiescent stage since the discov- 
ery observation. It is worth noting that none of the 
GC qLMXBs discovered in quiescence so far have been 
seen in outburst. Moreover, an outburst happening be- 
tween the available observations would have had an im- 
pact on the observed luminosit y and intensity variabil- 
ity w o uld have been d etecte d |Ushomirskv fc Rutledgd 
l200lt iRutledge et "aD I200I . The search of short- 
timescale variability (< 1 day) did not reveal any such 
variability. Finally, a PDS analysis was performed and 
did not demonstrate evidence of periodic variability in 
the frequency range 0.0001-0.10 Hz. The lack of inten- 
sity variability on various timescales further supports the 
classification of the source. 

In conclusion, the qLMXB U24 in the GC NGC 6397 
adds to the list of qLMXBs suitable to place constraints 
on the dense matter EoSs with a best-fit projected radius 
of i?oo = 11.9to^km. 
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